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Abstract 

We present measurements of decay matrix elements for hadronic transitions of the form T(nS) — > 
T(mS)-7nr, where (n, m) = (3, 1), (2, 1), (3, 2). We reconstruct charged and neutral pion modes with 
the final state Upsilon decaying to either [i + or e + e~. Dalitz plot distributions for the twelve 
decay modes are fit individually as well as jointly assuming isospin symmetry, thereby measuring 
the matrix elements of the decay amplitude. We observe and account for the anomaly previously 
noted in the di-pion invariant mass distribution for the T(3S) — > T(1S)tttt transition and obtain 
good descriptions of the dynamics of the decay using the most general decay amplitude allowed by 
partial conservation of the axial-vector current (PCAC) considerations. The fits further indicate 
that the T(2S) — ► T(1S)tttt and T(3S) — ► T(2S)mr transitions also show the presence of terms in 
the decay amplitude that were previously ignored, although at a relatively suppressed level. 

PACS numbers: 13.20.Gd,13.25.Gv,14.40.Gx 
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I. INTRODUCTION 



The transitions T(nS) — ► T(mS)7nT are of particular interest as probes of heavy quark 
and low energy QCD systems. The large b quark mass causes the bb bound state to have a 
very small radius (~ 1 GeV" 1 ) and to be non-relativistic ((v/c) 2 ~ 0.1). This makes these 
transitions ideal to study the process by which a pion pair is excited from the vacuum by the 
gluon field. The transitions among the massive bound states making up the T(nS) family 
can be calculated in terms of multipole moments of the chromo-dynamic field, providing 
simple relative rate and transition rule predictions. The pion pair excitation can be factored 
out and approximated separately. Most recent theoretical work has concentrated on this 
latter aspect of the decays. 

The f]' — > rj-KTi transition was the first decay of this form studied followed some 
years later by the if)' — *> J/ipirK transition 0. The rj' decay is only barely above threshold, 
and so the transition cannot show significant structure. Detailed study of the kinematics 
confirmed this. In contrast to this, the ip' decay has decay dynamics very different from 
a phase space distribution. The di-pion invariant mass distribution of this decay shows 
strong enhancement at larger values of M n7T . However, this is consistent with the presence 
of only the simplest term in the general Lorentz invariant amplitude derived from PCAC 
considerations (H, 0]. This is supported by the isotropic decay angular distribution of the 
pions, implying a minimal D-wave component. 

Previous CLEO data have been used to study T(nS) — > T(mS)7T7i transitions 0, 0, 0, 
with the T(2S) — > T(15')7r7r and T(3S') — > T(2S)tttt transitions following this same 
pattern in the di-pion invariant mass spectra as for the lighter mesons. But the T(3S) — > 
T(lS)inr transition has a second, strong rate enhancement near the tttt invariant mass 
threshold. This enhancement and the accompanying depletion at intermediate invariant 
mass are inconsistent with either pure phase space or the simple matrix element describing 
the ip' — > J I ipirn observations. Either another term must be included in the Lorentz invariant 
matrix element, or one must question the applicability of PCAC to the pion excitation and 
the validity of the multipole expansion of the bb bound state. 

Various mechanisms have been suggested to explain this anomaly, such as (i) large con- 
tributions from final state interactions 0, E3]j__(ii}a a isoscalar resonance in the mr sys- 
tem [ll|, (iii) exotic T — 7r resonances |9l [l3l 14 Jl5| . (iv) an ad hoc constant term in 
the amplitude [lg|, (v) coupled channel effects 1 7 ~l 1 8t| . (vi) S — D mixing [l9|, and (vii) 
relativistic corrections (20| . 

More recent experimental analyses with the very large data sets accumulated by the 



B factories at the T(4S*) show interesting behavior as well. Belle [2l| and B&B&r [22J do 
not see such anomalous behavior in the T(4S*) — > T(l S)ittt transition, but B&Bax does see 
such a double peaked structure in the T(4S') — > T(2S')7r7r transition. 

The shapes of the decay distributions originate in the details of the excitation of the pion 
pair from the vacuum and the particular projection of the initial state onto the final state. 
Hence, the enhancement of the decay rate at low M n7T , thus far considered an anomaly, is a 
good probe of the details of low energy QCD in the transitions of the bound states and the 
excitation of light hadrons from the vacuum. 

The general matrix element constrained by PCAC was derived by Brown and Cahn [3[ 
and is further constrained by treating the Upsilon transition as a multipole expansion as 



derived by Gottfried [231 ]. Yan [24J, Voloshin and Zakharov [25j, and others. The general 
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transition amplitude is then given in non-relativistic form: 



M = A(e' ■ e)(q 2 - 2M 2 ) + B(e' ■ e)E 1 E 2 + C((e' • 9l )(e • q 2 ) + (e' • g 2 )(e • q,)) , (1) 

where e' and e are the polarization vectors of the parent and final state Upsilons, and qi j2 are 
the four-momenta of the pions. In the first term, q 2 is the invariant mass of the pion pair. The 
quantities E\ and E 2 are the energies of the two pions in the parent rest frame, essentially 
indistinguishable from the lab frame due to the large masses of the Upsilons. 1 The third, 
or "C" term in this expression couples transitions via the chromo-magnetic moment of the 
bound state b quarks, hence requiring a spin flip. This is expected to be highly suppressed by 
the large mass of the b quark, so we expect only the first two terms to contribute. Neglecting 
the dependence on the parent and final state Upsilon polarizations (which apply only to the 
C-term), we have only two degrees of freedom, the Dalitz variables q 2 = M 2 n and r 2 = My w . 
In writing this amplitude, we have assumed the chiral limit, so that a fourth term, gM 2 , is 



taken to be zero |26|. 27 



The expression in Eqn. [T] can be made fully Lorentz invariant by rewriting the energy 
product in the B term as 

E\E 2 w l(P' ■ qi )(P ■ q 2 ) + (P' ■ q 2 )(P ■ gi )]/[2M T ,M T ], (2) 

with P' and P being the initial state and final state T four-momenta. 

The quantities A, B, and C are form factors that depend on the detailed dynamics of the 
decay. They are in principle functions of the Dalitz variables q 2 and r 2 . However, we expect 
them to vary on the scale of Aqcd, which is comparable to the total energy release of the 
decays, so to first order we assume they are complex constants. Angular structure or M nn 
dependence beyond that indicated in the explicit amplitude, Eqn. [U would be an indication 
of the non-constancy of these form factors, or alternately the breakdown of the assumptions 
leading to Eqn. Q3 

The di-pion transition can be interpreted as taking place in sequential two-body decays 
through a fictitious intermediate state X via the chain T(nS) — > T(mS)X and X — > Tin 
(see Fig. [p. In this view we can define the helicity angle of the X — > mi decay in the usual 
manner of the Jacob and Wick formalism. The polar helicity angle is referred to as Qx- 
Its cosine is used interchangeably with the second Dalitz variable, r 2 , as they are linearly 
related: 



2r 2 = 2M 2 + My, +M 2 -q 2 - cob8 x ^ (q 2 ~ 4M2) A 3 (M 2 ,, M 2 , q 2 ) (3) 

where A3 (a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 26c. These variables (r 2 or cos#x) carry 
structure from the second term in the amplitude due to the following relation: 

E 1 E 2 = \ ([E x + E 2 ) 2 - AE 2 max cos 2 6 X ) , (4) 



For transitions from the T(3S), the parent frame and lab frame are virtually identical. Even for T(2S) — ► 
T(lS , )7T7r transitions, in which the T(2S) comes from hadronic or electromagnetic transitions from the 
Y(3S), the parent's motion in the lab frame is unobservable other than in a small broadening of recoil 
mass peak and a slight smearing of reconstructed variables. 
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FIG. 1: (Left) The decay T(nS) — > Y(mS)irir follows the production of an initial state labeled 
T' which decays to an Yirir state. In our analysis, the final state T decays to a lepton pair whose 
momentum vectors are very nearly back-to-back due to the large energy release. (Center) The 
decay of the initial state Upsilon is governed by two kinematic variables, the Dalitz masses Mxn 
and M^. (Right) Alternately one can think of the irir system as a composite, X, and study its 
structure via the pion "decay" angles. 



with AE = E2 — E\. Because the initial state and final state Upsilons are essentially at 
rest, the energy sum E\ + E% is nearly a constant and equal to the mass difference between 
the Upsilons. For the 7r + 7r~ final state, Ox is defined as the angle of the positive pion, with 
— 1 < cos9 x < 1; for the 7r°7r° final state, because one cannot distinguish between the two 
neutral pions, we take < cos6>x < 1. 

Finding the presence of a non-zero C term would indicate the breakdown of the multi- 
pole expansion, i.e., of the assumption that the pion pair excitation is independent of the 
Upsilon transition process from n 3 Si state to m 3 Si, and that the spin flip of the b quarks is 
suppressed. However, finding a non-zero C term could also be due to distortions of the dis- 
tribution not accountable for by using only the first two terms with complex, but constant, 
coefficients A and B. 



II. DATA SETS AND EVENT SELECTION 



Data were collected with the CLEO III detector which is described in detail elsewhere [28 , 
29|, M- In this analysis we observe e ± , /j, ± , ir ± , and 7 particles in the final state, and 
so use both the tracking and calorimetry information from the detector, as well as lepton 
identification. Thus we employ global event, track, lepton, shower, and neutral pion selection 
criteria, in addition to signal and background identification criteria. 

The data were taken while running on the T(3>S) resonance, subject to standard CLEO 
data quality selections, and represent an integrated luminosity of 1.14 fb _1 , and an T(3S f ) 
production yield of (4.98 ± 0.01) x 10 6 . The T(2S) sample is obtained by reconstruction 
of sequential decays, T(35) — ► T(2S) + anything, occurring in this sample. The T(2S) 
population of (5.27±0.40) x 10 5 is estimated from the branching fraction [31] of 10.6%±0.8% 
for the decay T(35) — > T(2S) + anything, which is dominated by pion pair transitions and 
sequential photon decays through the Xfe(2P) states. 

All integrals needed in the analysis (for evaluation of acceptances and efficiencies) are 
calculated via the Monte Carlo method. Physics event generation is performed using the 



Lund Monte Carlo [32fl embedded in the CLEO physics Monte Carlo QQ [33j. The Lund 



event generator is used because it accurately accounts for the physics of the QCD bound 
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FIG. 2: Di-lepton invariant mass distributions for lepton pairs; the abscissa is the di-lepton 
invariant mass, showing peaks at the masses of the T(15) and T(25) mesons. The hatching 
indicates the limits to the invariant mass selection windows. Candidates are plotted after the 
signal selection described in Section II. D. At left are the di-muon candidates and at right the 
di-electron candidates. 



state production. The T(3S) produced in the e + e~ collision is then decayed according to 
standard decay tables and the detector response to the decay products is simulated using 
the physics simulation package GEANT 

In general, since all integrals are performed with respect to the natural measure over 
phase space, only phase space decays need be simulated. The decay amplitude is known 
exactly as a function of the decay kinematics, so all inputs to the matrix element extraction 
(other than acceptance and efficiency) are known to the precision of detector reconstruction. 

We select events containing two leptons or e + e _ ) and two pions {n + n~ or 7r°7r°). 

All low momentum tracks are assumed to be pions, because there is insufficient phase space 
for the production of a pair of kaons in a transition among any two of the three bound state 
Upsilons. Electrons and muons are identified by their energy loss and penetration depth in 
the detector as detailed below, and are required to be consistent with originating from either 
an T(2S) or an Y(IS') decay. The pion candidates are constrained to come from a common 
point at the beam location and the recoil mass (M r 2 ec = P rec • P rec ; P rec = Pb ea m — Qi — 12] 
see below) is used to identify the transition. The lepton pair invariant mass spectra and the 
recoil mass spectra are shown in Figs. [2] and [31 respectively. 



A. Global Event Selection 

The data used in this analysis are required to have been taken while running on the T(3S) 
resonance energy. Global event characteristics are used to preselect the events. Excessive 
tracks or showers in an event can dramatically increase the combinatoric background. To 
avoid this, reconstructed events are selected subject to upper limits on number of charged 
particle tracks and number of calorimeter showers. To establish conservative limits, signal 
Monte Carlo is studied for T(2S) — > T(1S)ttti transitions, which are the "worst case", in 
that extra tracks and showers in these modes arise from the initial transition from the T(3S) 
to the T(2S). Neglecting stray particles and secondary showers, there should be no more 
than four low momentum charged particle tracks and no more than eight electromagnetic 
showers in signal events. Comparison between data and Monte Carlo show good agreement 
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FIG. 3: Recoil mass, M rec , distributions for all modes. The upper plot is generated from neutral 
decays, T(nS) — ► T(mS)ir°ir° and the lower from charged decays, T(nS) — > T(mS)ir + ir~ . The 
final signal selections (track quality, pion quality, di-lepton mass, etc.) have been applied. The 
peaks at the T(1S) and T(2S) masses correspond to decays to these resonances from an T(3S) 
parent. The peaks at 9.8 GeV/c 2 are from T(2S)— > Y(lS)7T7r decays. The hatching shows the 
bounds on the recoil mass values for the three transitions. See also the window definitions in 
Tab. HI Yields are set to zero in the regions that correspond neither to signal nor to sidebands. 



in the number of tracks and showers found in the selected events. 



B. Selection of Final State Particles 

All candidate charged tracks are required to satisfy quality criteria. They must: 
-come from within 5cm of the origin along the beam axis (detector z axis); 
-come within 5mm of the beam axis (impact parameter); 
- have momentum less than the beam energy; 
-have a good helix track fit, with \ 2 P er hit less than 20. 
These requirements are applied to all track candidates and are augmented with identification 
criteria for leptons (see below) before being accepted as decay candidates. 

The charged transition pions frequently are of such low transverse momentum that they 
make two or more semi-circular arcs in the tracking volume. These "excess" tracks are 
removed by comparing the helix parameters, taking into account the expected energy loss 
as these pions spiral through the drift chamber. 

Candidate muons and electrons are required to have high momentum by requiring their 
transverse momentum to be px > 1 GeV/c, which removes a large fraction of the events 
with non-leptonic Upsilon decays. Because the leptons we seek originate from the decay of 
objects more massive than 9.4 GeV/c 2 , they pass this requirement easily. 

Muons are selected from among good tracks and are additionally required to penetrate 
the muon chambers to a depth of at least three interaction lengths. The ratio of energy 
deposition in the calorimeter to track momentum must also be less than one half, E/pc < 0.5. 

Electrons are selected from among good tracks and are additionally required to have a 
ratio of energy deposited in the calorimeter to track momentum E/pc > 0.5, as well as 
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having a profile of energy deposition consistent with that of an electromagnetic shower and 
a good spatial match between the shower and the track. The E/pc ratio selection is a very 
loose requirement added only as a precaution against muons contaminating the electron 
sample. 

The di-lepton mass is loosely required to be that of the final state Upsilon being studied, 
as shown in Fig. [2J For the T(1S) we require 9.25 < M« < 9.75 GeV/c 2 , while for the 
T(2S) we demand M u > 9.85 GeV/c 2 . Due to the large widths of these invariant mass 
peaks, no side band selection is performed in this variable, but rather only in the recoil mass 
distribution. 

The 7T° candidates are reconstructed from photon pairs. This begins by applying selection 
criteria to the showers. To be considered a photon, a shower must: 

- have energy greater than 30 MeV; 

- have a lateral shower profile consistent with that of a photon; 

- be inconsistent with the extrapolation of any track in the detector; 
-not include noisy channels in the calorimeter; 

- not be in the overlap region between the barrel and endcap calorimeter modules; 

- not be in the ring of crystals closest to the beam axis. 

Showers satisfying these selection criteria are considered to be photons and are combined 
into 7T° candidates. Photon pairs are required to have an invariant mass within 50MeV/c 2 
of the nominal 7r° mass, M n o. They are then required to fall within the asymmetric window 

_ 4< M„-M„ <3 (5) 

The photon-pair mass resolution, a 77 , is typically 5-7 MeV/c 2 . Candidate photon pairs 
are then kinematically constrained (subject to the measured uncertainties on energies and 
shower spatial locations) to have an invariant mass equal to M n o. To be used, 7r° candidates 
are further required to have a successful kinematic fit with confidence level (one degree of 
freedom) greater than 0.1%. 



C. Recoil Mass and Signal and Background Regions 

We select events for each transition by cutting on the mass of the system recoiling against 
the two pions in the T' — > 7r7r+" anything" decay: M 2 ec = M\, + q 2 — 2q ■ P', where, as 
above, q = q\ + q2 and P' is the Lorentz momentum of the initial state Upsilon. Given 
the large mass of the initial state Upsilon, the dot product simplifies and the recoil mass 
can be well approximated by M 2 ec « M\, + q 2 — 2Mt'(£>i + E 2 ). For the cascade decays, 
T(2S) — > T(lS')7r7r, this is not quite correct because the Lorentz momentum of the initial 
state Upsilon (the T(2S)) is not equal to the beam momentum. However, because the 
total momentum of the pions is small and the initial state is approximately at rest, using 
the incorrect momentum for the initial state does not significantly change the recoil mass 
distribution other than to shift it by the difference between T(3S) and T(2S) masses. Hence, 
we expect to find three recoil mass peaks. The transitions originating from the T(3S) 
will generate recoil mass (M rec ) peaks at the masses of the T(15) and T(2S), while the 
T(25) — > T(1S)titt decays will yield a peak at 9.79GeV/c 2 . These three peaks are clearly 
visible in Fig. [3j 
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FIG. 4: Candidate events that passed all selection criteria, and that have the final state Upsilon 
decaying to In the middle is the decay T(35) — ► T(lS)ir + ir~ . To the left is its neutral 

counterpart T(35) — ► T(1S)tt°tt . To the right is the charged transition T(25) — ► T(lS)ir + ir~ , 
with the bulk of its distribution at large values of dip ion invariant mass. In each plot there are ten 
degrees of grey-scale ranging from white (lowest occupancy per bin) to black (highest occupancy). 



The recoil mass, M rec , is measured rather accurately, especially in the charged case, due 
to the good resolution on the momenta of the low-momentum pions. It is still quite good 
for the neutral modes where the total pion momentum is given as the sum of momenta of 
two 7T° candidates reconstructed from the calorimeter showers. 



D. Signal and Background Selection 

The fit requires signal and background samples. They are determined as a function of 
the recoil mass, M rec , only. The recoil mass peak widths are determined from Monte Carlo 
with tight selections on variables other than the recoil mass. These widths are then used to 
determine mass windows to select events in both Monte Carlo and data samples. The signal 
regions are defined as the range within three times the peak width of the nominal recoil 
mass, while the backgrounds are the regions from six to twelve times the peak width from 
the nominal mass above and below the peak mass. The masses and widths used to define 
these regions are listed in Table [H The width of the recoil mass distribution in the decays 
T(2S) — > T(lS)n + n^ is roughly twice that of the direct decays. This is due to the boost 
of the initial state Upsilon imparted in its production by the cascade from the T(3S). The 
edges of the signal windows are indicated by the hatching in Fig. [3J Note that in Fig. [3] the 
yield in the regions not used for either signal or background definition have been set to zero. 

The Dalitz plot distributions for the selected data in six of the twelve final states are shown 
in Figs. H]and [5j Comparison of the 7r°7r° and 7i + tt~ for the T(35) — > T(1S)7ttt shows the 
depletion in charged particle efficiency at moderate di-pion invariant mass and large |cos 9x\- 
Comparison of the charged modes for T(35) — ► T(lS')7r + 7r _ and T(2S) — > T(lS)7i + 7i~ 
shows, in two dimensions, the obvious disparity between the two distributions. 
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FIG. 5: Candidate events that have passed all selection criteria, and that have the final state 
Upsilon decaying to e + e~. As in the prior plot, three transitions are, left to right, T(35) — > 
T(l 1 S)7r vr , T(35) T(1S)tt+it-, and T(25) -> T(lS)7r+7r-. 



Transition 


Recoil Mass Width (Data) Width (MC) Width (Cut) 






(MeV/c 2 ) 


(MeV/c 2 ) 


(MeV/c 2 ) 


(MeV/c 2 ) 


T(35) - 


-» T(15)vr+^- 


9 460.4 


2.4 


2.5 


2.5 


T(2S) - 


-» T(1S)tt + tt- 


9 792.4 


5.0 


5.0 


5.0 


T(3S) - 


-> T(25)tt+^- 


10 023.3 


2.2 


1.9 


2.1 


T(3S) - 


-» T(15)7r°7r° 


9 460.4 


15.0 


12.7 


13.8 


T(25) - 


-» T(15)vr vr 


9 792.4 


10.9 


10.5 


10.7 


T(3S) - 


-» T(25)7r°7r° 


10 023.3 


3.4 


3.4 


3.4 



TABLE I: Recoil mass distribution central values and widths for the signal and background 
selections used in the fit. The central values and the widths agree well between data and Monte 
Carlo. The signal windows are defined as the region within three times the cut width (last column) 
of the central mass and the background windows are defined as the region from six to twelve 
cut widths from the center on either side. The background subtraction is only important for the 
cascade decays for which there is a large contribution to the signal region from event combinatorics. 

III. MATRIX ELEMENT FITS 
A. The Likelihood Fitter 

The binned likelihood fit to the kinematic distributions of the T(mS)— > T(ns)7T7r decays 
is designed to deal correctly with the low bin yields expected from dividing approximately 
2000 events over a two dimensional space with more than ten bins per dimension. The 
general case of this problem is solved in Ref. [35]. Specific details of our application of this 
technique, including notes on variable smearing and background inclusion, are found in the 
Appendix. We fit the decay distributions to a product of the squared modulus of the decay 
amplitude and the phase space density sculpted by the detector acceptance. The matrix 
element has a known analytical form (see Eqn. [TJ as a function of the form factors A, B, 
and C, which are taken as complex constants. Its leading angular structure is known, and 
so long as the form factors are known, too, the entire amplitude can be described exactly. 
However, we cannot model the detector acceptance in analytic form, so we approximate its 
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FIG. 6: The efficiency-sculpted phase space in the two-dimensional plane for the transitions 
T(3S)-> T(lS)7r+7r- (left) and T(3S)-> T(lS)vr 7r° (right). Note that the neutral final state has a 
more uniform efficiency, especially in the region of moderate di-pion mass and large |cos(#x)|- 




FIG. 7: The three functions used in the fit for the T(3S) decay to T(1S)7T7t. From left to right 
these are for the pure A term, the interference term, and the pure B term. 

effect via Monte Carlo integration. 

We determine the integral of the phase space density in a bin in (g 2 , cos 9x), sculpted by 
acceptance and efficiency, by counting Monte Carlo events that pass the selection criteria and 
fall into that bin. In Fig. [6] we show the two-dimensional phase space after such sculpting. 
Note that while the overall efficiency for the neutral final state is lower than for its charged 
counterpart, the former is more uniform, particularly in the regions of intermediate M n7T and 
large |cos 9x\- For each bin of the observed distribution we predict the number of events 
as a function of the matrix element parameters by multiplying the Monte Carlo integral for 
that bin by the exactly calculated matrix element value for that bin. This approach avoids 
generating Monte Carlo integrated templates for each component of the angular distribution 
and reduces the uncertainty due to finite Monte Carlo sample size. 

To fit the decay distribution we take the squared modulus of the decay amplitude, Eqn. [[J 
and decompose it as a sum of six functional forms each multiplied by one of \A\ 2 , \B\ 2 , \C\ 2 , 
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$l(A*B), $l(A*C), or 3?(B*C). For normalization, the matrix element A is set to unity. 

The functional forms (e.g., (q 2 — 2M 2 ) 2 ) depend on the Dalitz variables and are pre- 
evaluated into templates over the Dalitz space. The fitter then seeks the best fit as a 
function of the matrix element ratios A, B, and C. The input to the fitter consists of only 
the data, background, and phase space Monte Carlo binned across the Dalitz plot, and the 
component templates of the decay distribution derived from the exact decay amplitude, but 
taking into account the kinematic smearing and acceptance and efficiency effects due to 
reconstruction as determined from the detector simulation. The backgorund component is 
scaled by the ratio of the signal region width (6 a; see Section II. D) to the total backgorund 
sideband width (nominally 12a). 

In Fig. [7] we show the functional forms for \ A\ 2 , $t,(A*B), and \B\ 2 for the case of T(3S)— > 
T(lS)7T7r. In our experiment, the complementarity of the neutral and charged final states is 
particularly important in that the rightmost of these (the form for |£>| 2 )) depletes the region 
for which the 7i + ir~ channel has falling efficiency. Consistent results between the 7r°7r° and 
7r + 7r~ transitions gives us confidence that the simulation of this fall-off in efficiency is reliable. 
The matrix element extraction procedure is tested "end-to-end" by simulating signal with 
known matrix elements in Monte Carlo and comparing the fit result and its uncertainty 
with the known inputs. Samples of the same size as the observed yield are generated and 
fit identically to the data. The results yield standard normal distributions in the observed 
uncertainty scaled residuals for widely distributed seed matrix element values. This confirms 
the fitter is unbiased at the level of precision to be expected from the sample size of the 
measurement. 

B. Fits with C = 

The fits to the two dimensional distributions of M n7T and cos^x determine the matrix 
element B/A and C/A. The extracted values of ^{B/A) and ^s{B/A) are summarized in 
Table [TT1 subject to the constraint that C = 0. In that we only measure the cosine of the 
phase difference between B and A, Q(B/A) is only known to within a sign. The upper set 
of matrix elements are obtained from independent fits to ten individual decay modes; we 
cannot individually fit the two modes associated with T(3S)^ Y(2S)7r + 7r~ because of their 
limited statistics. The lower set of three are from the simultaneous fits of all final states for 
each given Upsilon transition. 

In the simultaneous fits the relative branching ratios between modes are not constrained, 
but it is assumed that the di-pion excitation dynamics is independent of the charge of the 
pion final state (isospin symmetry) and thus the decay distributions should be identical 
to within statistical fluctuations for all transitions between the same Upsilon states. This 
assumption is supported by the consistency among the matrix element values extracted 
independently, as well as their consistency with the value extracted from the simultaneous 
fit. In particular, the four final states studied for the transition from T(3S) to T(1S) show 
excellent agreement between the two lepton species and between charged and neutral pions. 

To study the fit quality we project the data and the expected decay distribution for the 
matrix element value preferred by the fit onto the di-pion mass (M n7r ) and di-pion helicity 
angle (cos#x) variables and calculate a x 2 f° r each projection. To increase the bin contents 
we sum over lepton species but not over pion charges. We expect the shapes for charged 
and neutral pions to differ due to the rather different efficiencies for reconstruction and 
resolutions, as well as the folding of the neutral angle in the fits. Figure M presents plots of 
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Individual Fits 








R(jB/A) 


T(35) 


-► T(15)vr+^- 


T- 




-2.514 ± 0.037 ±1.164 ± 0.059 






T - 






-2.527 ± 0.049 ±1.180 ± 0.079 


T(35) 


T(15)7r°^°; 


T - 






-2.426 ± 0.085 ±1.313 ± 0.159 






T - 






-2.524 ± 0.093 ±1.070 ± 0.153 


T(25) 


-► T(15)tt+^- 


T- 




V- 


-0.656 ± 0.126 ±0.431 ± 0.089 






T - 




~e~ 


-0.689 ± 0.147 ±0.425 ± 0.102 


T(25) 


-» T(15)7r°7r°; 


T - 


-> /I' 


V" 


-0.148 ±0.280 0.000 ±1.655 






T - 
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~e~ 


-0.293 ± 0.330 ±0.001 ± 1.130 


T(35) 


-» T(2S')7r 7r ; 


T - 


-» M" 


v- 


-0.283 ± 0.305 ±0.001 ± 1.708 
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Simult 


aneous Fits 








U(B/A) $S(B/A) 


T(35) 


-» T(15)vr^ 








-2.523 ± 0.031 ±1.189 ± 0.051 


T(25) 


-» T(15)vr^ 








-0.753 ± 0.064 ±0.000 ± 0.108 


T(35) 


-» T(25)vr^ 








-0.395 ± 0.295 ±0.001 ± 1.053 



TABLE II: Fit results from T(nS) — ► T(mS)irir transitions for B/A with C set to zero. The upper 
set of results is from individual fits to each separate decay mode and the lower set of results is 
from simultaneous fits to both lepton final states and both pion charge modes. We cannot fit the 
Y(35) — > Y(25)"7r + 7r~ transitions, individually in e + e~ and fi + f^ or combined, because of their 
limited statistics. In the simultaneous fits the relative branching fractions are allowed to float. 
Note that we know the value of the imaginary part of the ratio only to within a sign. 



Upsilon 
Transition 


7T + 7T 


ttV 


cos Ox 


M W7T 


cos 6x 


M W7T 


35 -» 15 
25 -» 15 
35 -+25 


33.2 (16) 
6.1 (10) 
7.1 (7) 


46.9 (32) 
22.7 (12) 
7.8 (6) 


4.3 (8) 

3.4 (5) 
7.4 (4) 


52.1 (32) 
13.7 (12) 
2.5 (7) ' 



TABLE III: The figure of merit for each of the twelve projections in the accompanying figure. For 
each projection we give the value of \ 2 and, in parentheses, the number of bins used to calculate 
it. Uncertainties in the fit results due to limited simulation statistics are not included in these 
calculations. 

the data overlaid with the fit results, showing good qualitative agreement. The \ 2 values 
from these overlays, given in Table IHIl are acceptable, given the simplicity of the fitted 
matrix element. 

As a further fit quality test, we examine the two dimensional distribution over the Dalitz 
variables of error- normalized deviations. The deviations, Si, are the difference, fit subtracted 
from the data, divided by the mutual uncertainty: 

5* = (6) 

where each di is the predicted decay population in bin i. The bin-by-bin uncertainties, Oj, 
are composed of the uncertainty on the data yield in the bin, = y/d~i, and the uncertainty 
on the template function, dominated by the fluctuation in the Monte Carlo phase space yield 
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FIG. 8: Plots overlaying projections of the data (points with error bars) and the fit result (his- 
tograms) onto the M n7T and cos Ox variables. The plots are summed over electrons and muons, but 
are differentiated by pion charge. The neutral modes (open symbols, dashed lines) show only a 
positive distribution in cos Ox because the two pions are indistinguishable. For the charged modes 
(solid symbols, solid lines) the angle is that of the tt + . 



and proportional to 1/y/ai, where dj is the Monte Carlo phase space yield in bin i. Hence, 
Oi = \Jdi + df/cii. 

The bins for which di — require special treatment, and cij is modified appropriately. To 
minimize the effect of such bins with zero yield, we sum over muon and electron final states. 
This takes a weighted average over the distributions, rather than taking account of the 
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M„ (GeV/c 2 ) 

FIG. 9: Plots of the bin by bin deviations of the data from the fit templates normalized to the 
expected uncertainty on the bin content for the transitions between T(3S) and T(1S). The left 
plot is for the charged pion modes while the right plot is for the neutral pion modes. The data are 
summed over lepton species. No strong concentration of deviations is apparent. 

differences between the individual distributions and their individual template predictions. 

The deviations between the data and the fit templates, Si, are shown in Fig. [9] for the 
charged and neutral transitions between T(3S) and T(1S). No significant bunching is ob- 
served that would indicate a bias. We neglect the small accumulations in the areas of low 
tracking efficiency (at large Icos^xl and intermediate M nn ), probably attributable to the 
Monte Carlo detector model not being sufficiently accurate. 

C. Fits Including the Chromo- magnetic Term C 

The fit results in Table [III do not take into account the possible presence of amplitude 
terms that come from chromo-magnetic couplings, which would allow the additional C term 
to appear. This term is nearly degenerate with the B term, and fits allowing it to float show 
a strong covariance between these two terms. This is caused by the similarity in structure 
of the two terms; B accompanies a functional dependence E X E 2 , while (e' • 91,2) (e • 92,1) 
emphasizes the regions of phase space in which the pion spatial momentum, and hence also 
the energy, are large. The low yield modes do not allow the measurement of the term at 
all. We therefore only study it in the T(3S) — ► T(lS)7nr transitions, and then only extract 
a value from the simultaneous fit. 

The covariance between B and C for the T(3S) — ► T(lS')7r7r transition is summarized 
in Fig. [TDJ which shows the variation of extracted |£>/.A| with |C/^4|, both as a fit error 
ellipse, and as fit trials with \C\ constrained to different values. The ellipse corresponding 
to one standard deviation from the best fit gives a value for T(35) — > T(1S)7ttt of \C/A\ = 
0.45 ± 0.18, with the uncertainty being purely the statistics of the fit. The fit which includes 
real and imaginary parts of C/A shows an improvement over the one with C fixed at zero of 
— 21n£ = 9.4. Although this implies a ~ 3a improvement in fit quality when C is allowed to 
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FIG. 10: Variation of B with C magnitudes. The points indicate the fit and error for B at fixed 
values of C. The ellipse indicates the one sigma bound on the free fit, the axis of which agrees well 
with the point by point fits. The bands indicate the one standard error limits on B when C is fixed 
to zero. 



float, systematic uncertainties, which are significant, have not yet been taken into account. 

With this extended fit the six projections of Fig. [8] show no significant changes, and for 
the T(3S) — > T(l S)ttit transition the best fit value of changes minimally from 2.79 

(C = 0) to 2.89 (C floating). The phase of B with respect to A, denoted 5ba, changes little 
(about 2 degrees) from the 155 degrees of the fit done with C — 0. The smallness of the 
effects is not surprising as the shapes of the B and C components of the amplitude are nearly 
degenerate. A non-zero value of \C/A\ may be a consequence of statistical fluctuations and 
small systematic biases or may be due to A and B having some dependence on q 2 and/or 



r , i.e. 



not being complex constants. 



D. Partial Wave Decomposition 

Since the focus of this study is the decay dynamics of the di-pion system it is useful to 
think about the spin structure of the di-pion composite. The idea is to look for signatures 
of higher spin resonances in the form factors A and B. We must account for the intrinsic 
spin structure of the Lorentz amplitude to do this. We equate the Lorentz amplitude with 
the general partial wave amplitude to relate the matrix elements. 

The transition is of the form (T; X\T'). If the di-pion system has spin J we have: 

(l,m r ;Jx,mx\l,m T ,) . (7) 
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In that here we assume that only A and B are non-zero, there is no change in the polarization 
from the initial state to final state Upsilon; more general partial wave decompositions can 



also be made [10l. 127]. The angular momentum projections are then mp = rat, and mx = 0. 
Hence the partial wave decomposition of the X system can only have m = components. 
Since the pions are in an iso-singlet state, their parities require their relative orbital angular 
momentum to be even, and hence the orbital angular momentum between the final state 
upsilon and the di-pion composite must also be even. We can only have even partial waves 
in our decomposition: 

M P = S(q 2 )Y°+V(q 2 )Y° 

= 5(g 2 )^ + P(g 2 ) V / J(|cos^ x -l). W 

The functions S(q 2 ) and T>(q 2 ) are composed of two terms each, one from the A dependence 
and one from the B dependence: 

S(q 2 ) = AS A (q 2 ) + BS B (q 2 ), and V(q 2 ) = AV A (q 2 ) + BV B (q 2 ). (9) 

We here assume that there are no significant contributions from partial waves higher than 
J = 2. This will be true if there are no contributions from variations of form factors over 
the Dalitz space. Higher J terms must originate from structure in the form factors A and 
B. 

Equating the decay distributions (or equivalently, projecting inner products over the 
angular space) yields the following forms: 

S A (q 2 ) = q 2 -2M 2 , and V A (q 2 ) = (10) 

for a pure "A" decay, and 




for a pure "£>" decay. The overall amplitude is 

M P = {AS A {q 2 ) + BS B (q 2 ))Y ° + (AV A (q 2 ) + BV B (q 2 ))Y 2 ° , (12) 

where it is implied that Y™ is a function of the helicity angles of the pseudo-decay X — > 7r7r, 
9x and <fix (although the latter variable plays no role in the description of this decay, by 
the assumptions above). Interference between the S- wave and D-wave components of the 
decay comes from the functions S(q 2 ) and T>(q 2 ) being complex valued. Though S A)B (q 2 ) 
and D AB (q 2 ) are real functions, A and B are complex coefficients with nontrivial relative 
phase. 

The structure of S and D components as functions of q 2 are determined by the assump- 
tions underlying the derivation of the general Lorentz amplitude. The four functions from 
the pure A and pure B components are sketched in Fig.[TT]together with the fractional S- and 
.D-wave components in the angular distribution (which can alternately be thought of as the 
strengths of the S- and D-wave components), extracted from our fit to T(35) — > T(lS)mi. 
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FIG. 11: The left plot shows the amplitude component functions 5^4, P4, and T>q as a function 
of = yfq 2 . These are summed to obtain the total amplitude. The partial rate to S-wave and 
D-wave components are shown in the right plot for the T(35) — > T(1S)ttit decay as determined 
from the results of this analysis: B/A = — 2.52+1. 19£. Note that the D-wave contribution is largest 
in the low to intermediate range of q 2 , and is suppressed at both extrema by angular momentum 
barrier effects. Note further that this is not a resonance phenomenon despite its shape in M n7T and 
the changing angular structure. 



This partial wave extraction becomes much more complex if the form factors are assumed 
to be variable over the Dalitz space, for example due to resonant structure/enhancement in 
the decay. This will introduce higher powers of cos 2 9x to the overall amplitude and will 
need higher partial wave components to account for the variation. 

The presence of D-wave components in the angular distribution of the decay is not in 
itself an indication of resonances contributing, nor the presence of unaccounted-for physics. 
The presence of a g 2 -dependent D-wave component could simply be a consequence of angu- 
lar momentum barriers in the three body phase space of the decay. The data do not demand 
the introduction of a g 2 -dependent magnitude or phase for B/B. These small D-wave com- 



ponents are consistent with those derived in a recent paper by Voloshin [20(], in which he 
emphasizes the importance of relativistic and chromo- magnetic effects. 



IV. SYSTEMATIC UNCERTAINTIES 

We address three sources of systematic uncertainty in the measurements of B/A and 
C/A: model dependence, detector efficiency and resolution, and backgrounds. 

In Sect. Ill we showed that our model provides a very good description of the data in the 
(q 2 , cos8x) plane and that the presence or absence of the chromo-magnetic coupled term in 
the amplitude has little effect on \B/A\ and Sba- 

Uncertainty in the estimation of the detector efficiency and resolution contributes most 
significantly in the charged mode analyses due to our limited knowledge of the tracking 
efficiency at very low momentum. In that the low momentum region is precisely where 
the matrix element has potential suppression in the B term, this can potentially cause a 
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significant bias. To estimate this effect we use the full Monte Carlo simulation with looser 
and tighter track reconstruction requirements to provide bounds on the shape of the efficiency 
as a function of track curvature. We then create a number of analytic functions that span 
these boundaries. Then we use a toy Monte Carlo to simulate events with one of these 
analytic functions and assume a different one for the reconstruction. The variations in the 
fit results are conservatively assumed to be one standard error uncertainties on the extracted 
parameters. 

The same process is repeated for the neutral modes, varying the thresholds at which show- 
ers can be observed in the detector. This obviously leads to a large variation in branching 
ratios from simple inability to reconstruct the decays, but does not exhibit any significant 
change in the shape of the efficiency function over the measurement variables. This is to be 
expected since the ir° decays have largely flat acceptance over the kinematic range of these 
decay modes. 

We have evaluated the systematic errors associated with detector resolution, and find 
them to be negligible in comparison with the statistical errors from the fit and the other 
systematic errors discussed here. The curvatures of the matrix element components across 
the Dalitz plot are all very much smaller than the variances of the reconstructed measurement 
variables around their true values. No systematic uncertainty is assigned to this source. 

Background subtraction is only a source of bias if the upper and lower sidebands in the 
recoil mass exhibit markedly different shapes or the background is strongly peaked under 
the signal. In this case the extrapolations of the background shape and magnitude under 
the peak could be distorted. We have redone the fits with the ratio of the widths of sideband 
window to signal window both doubled and halved, and with only using either the high-mass 
or low-mass sideband. The variations in the fit are conservatively taken to represent one 
sigma variations in the final result, and are given in the last column of Table HVl 

Finally, the lepton reconstruction is capable of contributing bias since all decay modes are 
fully reconstructed. However, the detector response to leptons is sufficiently well measured 
in other analyses that the detector simulation is much more precise than what is required 
for this data set. The variation of the shapes is furthermore only relevant for the final C 
term, which is dependent on the lepton polar angle. With the exception of a small part of 
the C terms there can be no effect due to lepton acceptance. We estimate any systematic 
error associated with the lepton reconstruction to be negligible. 

The fit results combined with these systematic uncertainties are summarized in Tables HVl 
and |Vl Since the magnitude \C/A\ in the fit is only separated from zero by about one 
standard error and is expected to be suppressed in the theoretical models, we set a limit 
rather than claim observation of a non-zero value. 

We set this limit by assuming the value of C/A has a Gaussian uncertainty in real and 
imaginary parts. We transform variables to and arg(C/^4), using the sum of the 

variances of statistical and systematic origin as the overall variance. We then find the 90% 
upper limit from the resulting distribution as 

\C/A\ < 1.09 at 90% C.L. . (13) 

V. SUMMARY AND ACKNOWLEDGMENTS 

We quote fit results for the three transitions from simultaneous fits to the different decay 
modes with statistical and systematic uncertainties in Table [V] Only the simplest features 
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Fit, No C 








stat. 


effcy. (tt^) 


effcy. (n° 


) bg. sub. 


T(35) -» 


T(15)vr^ 


3(B/i) 


-2.523 
±1.189 


±0.031 
±0.051 


±0.019 
±0.026 


±0.011 
±0.018 


±0.001 
±0.015 


T(25) -» 


T(15)vr^ 


U(B/A) 
%{B/A) 


-0.753 
0.000 


±0.064 
±0.108 


±0.059 
±0.036 


±0.035 
±0.012 


±0.112 
±0.001 


T(35) -» 


T(25)vr^ 


SO(U 1 A\ 

S(B/A) 


— u.oyo 
±0.001 


_i_n oq^ 
±1.053 




_i_n no^ 
±0.180 


_i_n i on 

ztU. 1ZU 

±0.001 


Fit, float C 






stat. 


effcy. (tt^) 


effcy. (vr° 


) bg. sub. 


T(35) -» 


T(15)vr^ 


\B/A\ 
\C/A\ 


2.89 
0.45 


±0.11 
±0.18 


±0.19 
±0.28 


±0.11 
±0.20 


±0.027 
±0.093 



TABLE IV: Combined fit results for all transitions with statistical and systematic uncertainties. 
The systematic uncertainties are in order: -k^ detection efficiency, tt° detection efficiency, and 
background subtraction for the T(3<S) — > T(2S)irir transition. The upper set of results are for 
the fits assuming contributions to the amplitude from only the A and B terms. The bottom two 
lines are the fit results when the C term is allowed to be non-zero. The imaginary part of the 
ratio has a two-fold ambiguity and is only known to within a sign. Note that for the transition 
T(3S) — > T(2S , )7T7r we do not have fits for the charged di-pion case. 



of the Brown and Calm decay amplitude (Eqn. [T]) are included in our model, and the fits 
account for the structure of the decay without introduction of new physics or contributions 
from resonances. 

The matrix elements are indicated as points in the complex plane in Fig. [121 For the 
"anomalous" T(35) — ► T(lS)mr transition we fit for the presence of the "suppressed" C 
term as a test for the breakdown of the underlying assumptions leading to the standard 
matrix element. This term is not significant when systematic errors are taken into account 
and the quality of the fit to the data is good without it. Therefore, we set an upper limit of 
\C/A\ < 1.09 at 90% C.L.. 

We note in particular that the treatment of the di-pion transitions via the full allowed 
matrix element under the assumptions in Refs. [1,0, 23, 24, 25] allows two matrix elements, 



only one of which has traditionally been assumed to be non-zero. The description of the 
T(3S) — > T(1S)tt7t transition di-pion mass and angular structure as anomalous is only true 
in the limit of this assumption. This analysis shows in particular that the description of the 
decay process in terms of the two favored amplitude terms, with complex form factors con- 
stant over the Dalitz plane, suffices to describe the decay distributions of T(3S) — > T(1S)tt7t, 
T(3S) — > T(2S)inT, and T(2S) — > T(1S)tttt, provided the form factors are allowed to vary 
with the transition. For the T(3S)^ T(1S)7T7t transition, we find |i3/^4| = 2.79±0.05, which 
could imply a large magnitude of B or a suppressed A; recent theoretical considerations [2p| 
favor the latter interpretation. While smaller than in the case of T(3S) — > T(lS)7nr, \B/A\ 
is also determined to be non-zero for the case of T(25') — > T(lS')7r7r. The large imaginary 
part of B/A is intriguing [27| . 

While there are not yet first principles predictions of the values of the matrix elements 
of the decays studied here, this analysis does provide complete measurements of the relative 
matrix element magnitudes and phases that can serve as a point of comparison with ab initio 
QCD calculations. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
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Fit, no C, total error 



T(35) -► T(lS)7T7r 



%1{B/A) 
%B/A) 

\B/A\ 

5ba 



-2.52 ± 0.04 

±1.19 ± 0.06 

2.79 ± 0.05 

155(205) ± 2 



T(25) -» T(15)7T7r 



I^MI 



-0.75 ± 0.15 

0.00 ± 0.11 

0.75 ± 0.15 

180 ± 9 



T(35) -> T(25)7T7r 



-0.40 ± 0.32 
0.00 ± 1.1 



Fit, float C, total error 



T(35) -» T(15)7T7r 



I^MI 
|CM| 



2.89 ± 0.25 
0.45 ± 0.40 



TABLE V: Fit results for all transitions with total uncertainties. These numbers represent the final 
result of this analysis. In the case of the magnitude ratio |C/«4|, we also quote a limit as detailed 
in the text. The phase angles are quoted in degrees, and have a two-fold ambiguity of reflection in 
the real axis. 
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FIG. 12: Complex values of matrix element ratio B/A from combined fits for the three transitions 
under the assumption that C = 0. Note the two-fold ambiguity in the imaginary part. 
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VI. APPENDIX: DETAILS OF THE LIKELIHOOD FITTER 



This appendix gives some details of our application of the likelihood fitter. 

Smearing due to reconstruction resolution adds a small variance to the Poisson error on 
the Monte Carlo integral, but the smearing widths are small compared to the scales over 
which the matrix element changes so this additional variance is small. For any shape with 
an approximately polynomial form at a point, the resolution is described by convolving a 
Gaussian with the polynomial. As an example, we assume a functional form g T = a+bx+cx 2 
and seek its observed shape in terms of the observed variables, g°(x°), using a Gaussian 
transformation: 

= J dx T g{x T -x°\n = 0, a)g T (x T ) (14) 

= J dx T g(x T -x°\fx = 0,a)(a + bx T + c(x T ) 2 ) (15) 
= (a + ca 2 ) + bx° + c(x°) 2 (16) 

So long as a 2 <C a/c, i.e., the resolution is small compared to the curvature, the shape 
will not be materially changed. For the angular dependence, which is quartic in cos 9x this 
means the resolution need only be small compared to 1/2; the observed resolutions are of 
the order of 5% or less. In M nn the same holds true, with the scale being given by the pion 
mass, 140MeV/c 2 , and the observed resolutions being at worst 10MeV/c 2 . The shape of 
the decay amplitude is not changed significantly by these resolutions, but any residual effect 
is included in the estimated tracking and shower systematic uncertainties. 

Our problem differs from that discussed in Ref. j35| in that the templates do not have 
independent Poisson fluctuations. The underlying phase space simulation has a Poisson 
fluctuation, but the templates are known (very nearly) exactly and uncertainties on them 
do not contribute to the overall likelihood function. 

In the absence of background this problem is solved as follows, with each two-dimensional 
(q 2 , cos6*x) bin denoted by subscript i. 

We compare the Monte Carlo simulated, acceptance and efficiency-corrected, phase space 
distribution (with true and observed yields A4 and a,), multiplied by the modulus squared of 
the amplitude, with the data distribution (with true and observed yields and dj). Both 
distributions are subject to Poisson fluctuation: 

V(d l ;D l )= 1 and P(a l ;A l ) = - (17) 

dil ail 

Bin-by-bin, the modulus squared of the decay amplitude appears in the exact relation be- 
tween the true data yields Di and the true phase space yields 

Di = fi{a)Ai. (18) 

The function represents the decay distribution (|.M| 2 ) in the kinematic space bin % as a 
function of a, the decay parameters. In this case a consists of real and imaginary parts of 
B/AandC/A. 

The log likelihood used in this fit is then given by, summing over all the bins: 

n 

\nC(a) = {dilnfi{a)Ai - fi{a)Ai - lndj + a t In A - A % - ma*!) . (19) 
i=i 
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The Ai represent the phase space subject to efficiency and acceptance effects and are 
uninteresting nuisance parameters that can be eliminated by extremizing the likelihood 



with respect to them. Proceeding in analogy with the approach in [35| we can find the 
analytic extremum condition, solve for Ai 

Ai = ^TT ( 2 °) 

Ji ~r J- 

and substitute back into the likelihood function to give a reduced likelihood: 

n r -i 

ln£(q) — di In fi(a>) — (di + aA ln(l + /»(«)) + const. (21) 

i=l L 

We then minimize —2 In £ with respect to the fit parameters a (occurring only in the co- 
efficients fi). This is implemented using the CERN Library minimization package, MI- 
NUIT HJ. 

The full likelihood as used in the fit includes an extension of this approach to account 
for background under the signal peaks. This introduces additional parameters Bi and 6,. 
These represent bin by bin true and observed background yields. The are a second set 
of nuisance parameters that are eliminated in the same way as were the Ai before. The 
resulting likelihood is significantly more complicated in detail but not in principle. For 
brevity it is not included here. 



23 



[1] A. Rittenberg, UCRL-18863, Ph.D. Thesis (unpublished), 1969; G. W. London et al, Phys. 

Rev. 143, 1034 (1966); J. Badier et al, Phys. Lett. 17, 337 (1965). 

[2] T. M. Himel, Ph.D. Thesis, SLAC-223 (1979); M. Oreglia et al. (Crystal Ball Collaboration), 

Phys. Rev. Lett. 45, 959 (1980); H. Albrecht et al. (ARGUS Collaboration), Zeit. fur Phys. 

C35, 283 (1987); J. Z. Bai et al. (BES Collaboration), Phys. Rev. D62, 032002 (2000). 

[3] L. S. Brown and R. N. Cahn, Phys. Rev. Lett. 35, 1 (1975). 

[4] M. B. Voloshin, P. Zh. Eksp. Fiz. 21, 733 (1975) [JETP Lett., 21, 347 (1975)]. 

[5] F. Butler et al. (CLEO Collaboration), Phys. Rev. D49, 40 (1994). 

[6] S. Glenn et al. (CLEO Collaboration), Phys. Rev. D59, 052003 (1999). 

[7] J. P. Alexander et al. (CLEO Collaboration), Phys. Rev. D58, 052004 (1998). 

[8] I. C. Brock et al. (CLEO Collaboration), Phys. Rev. D43, 1448 (1991). 

[9] G. Belanger, T. DeGrand and P. Moxhay, Phys. Rev. D39, 257 (1989). 

[10] S. Chakravarty, S. M. Kim and P. Ko, Phys. Rev. D50, 389 (1994). 

[11] T. Komada, S. Ishida and M. Ishida, Phys. Lett. B 508, 31 (2001); Phys. Lett. B 518, 47 
(2001). 

[12] M. Uehara, Prog. Theor. Phys. 109, 265 (2003). 

[13] M. B. Voloshin, P. Zh. Eksp. Fiz. 37, 58 (1983) [JETP Lett. 37, 69 (1983)]. 

[14] V. V. Anisovich, D. V. Bugg, A. V. Sarantsev and B. S. Zhou, Phys. Rev. D51, 4619 (1995). 

[15] F.-K. Guo, P.-N. Shen, H.-C. Chiang and R.-G. Ping, Nucl. Phys. A 761, 269 (2005). 

[16] P. Moxhay, Phys. Rev. D39, 3497 (1989). 

[17] H. J. Lipkin and S. F.Tuan, Phys. Lett. B 206, 349 (1988). 

[18] H.-Y. Zhou and Y-P. Kuang, Phys. Rev. D44, 756 (1991). 

[19] S. Chakravarty, S. M. Kim and P. Ko, Phys. Rev. D48, 1212 (1993). 

[20] M. B. Voloshin, Phys. Rev. D74, 054022 (2006). 

[21] K. Abe et al. (BELLE Collaboration) (2006), hep-ex/0611026. 

[22] B. Aubert et al. (B&Bm Collaboration), Phys. Rev. Lett. 96, 232001 (2006). 

[23] K. Gottfried, Phys. Rev. Lett. 40, 598 (1978). 

[24] T. M. Yan, Phys. Rev. D22, 1652 (1980). 

[25] M. B. Voloshin and V. I. Zakharov, Phys. Rev. Lett. 45, 688 (1980). 

[26] T. Mannel and R. Urech, Z. Phys. C73, 541 (1997). 

[27] S. Dubynskiy and M. B. Voloshin. la7Xiv:0707. 12721 

[28] G. Viehhauser, Nucl. Instrum. Meth. A462, 146 (2001). 

[29] D. Peterson et al, Nucl. Instrum. Meth. A478, 142 (2002). 

[30] M. Artuso et al, Nucl. Instrum. Meth. A554, 147 (2005). 

[31] W.-M. Yao et al. (Particle Data Group), Jour, of Physics G 33, 1 (2006). 

[32] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 (2001); for more details see T. Sjos- 

trand, L. Lonnblad and S. Mrenna, PYTHIA 6.2 Physics and Manual, |hep-ph/01 08264, 

[33] See http : //w4 . Ins . Cornell . edu/public/ CLEO/sof t/QQ/index . html, 

[34] GEANT: Detector Description and Simulation Tool, CERN (1993), long writeup W5013. 

[35] R. J. Barlow and C. Beeston, Comput. Phys. Commun. 77, 219 (1993). 

[36] F. James, MINUIT: Function Minimization and Error Analysis, CERN (1998), long writeup 
D506. 



24 



